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Abstract 

The interactions of Z)S C -Z)A C , D*T, C -D* A c , and related strangeness channels, are studied within 
the framework of the coupled channel unitary approach with the local hidden gauge formalism. A 
series of meson-baryon dynamically generated relatively narrow A'"* and A* resonances are predicted 
around 4.3 GeV in the hidden charm sector. We make estimates of production cross sections of 
these predicted resonances in pp collisions for PANDA at the forthcoming FAIR facility. 
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I. INTRODUCTION 



The use of chiral Lagrangians in combination with unitary techniques in coupled chan- 
nels of mesons and baryons has been a very fruitful scheme to study the nature of many 
hadron resonances. The poles found in the analysis of meson baryon scattering amplitudes 
are identified with existing baryon resonances. In this way the interaction of the octet of 
pseudoscalar mesons with the octet of stable baryons has lead to J/P = l/2~ resonances 
which fit quite well the spectrum of the known low lying resonances with these quantum 
number [lH5|. The combination of pseudoscalars with the decuplet of baryons has also 



received attention and also leads to several dynamically generated states 



Work sub- 



stituting pseudoscalar mesons with vector mesons has also been done recently leading to 



new resonances dynamically generated 



9]. 



One of the interesting findings in the study of the interaction of pseudoscalars with the 
octet of baryons is the generation of the iY"*(1535) resonance which has large couplings to 
and KA, to the point that the resonance can be approximately considered as a bound state 



of these meson baryon components 



12Hl4j . Another point of view is that this resonance can 



be considered as a hidden strangeness state. In fact, phenomeno logical studies show that, 
indeed, this seems to be the case 



15 



16|. 



The idea that we want to explore here is to see if one can also generate dynamically 
baryon states in the hidden charm sector. The interaction of charmed mesons with the octet 



of stable baryons has been studied in 17|, |l8| and further refined in 19h21| . Several states 



with open charm are dynamically generated there, in particular the A c (2593). 

In the present work we follow the steps of flQent eeneentrate in states ef hidden ehann, 
for which we study the interaction of an anticharmed meson with a charmed baryon. The 
und erly ing; theory that we use is an extension to SU(4) of the local hidden gauge Lagrangians 



22h25|. where SU(4) is broken to account for the different masses of the vector mesons 



exchanged in the t- and u- channels. The study is done both with pseudoscalar mesons and 
vector mesons and we obtain three dynamically generated hidden charm baryons generated 
from the pseudoscalar baryon interaction plus three other states from the interaction of 
vector mesons with baryons, all of them with masses around 4200-4600 MeV. 

We also make estimates of production cross sections with p collisions that could be carried 
out at the future FAIR facility within the PANDA project. We also study how the presence 
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of these resonances could increase the rate of J/ip and r\ c production around the energies 



26|, 



where the resonances can be formed. Part of our results have been briefly reported in 
here we give a much more complete report of our investigation. 

In the next section, we present the formalism and ingredients for the study of the in- 
teraction, and give the poles obtained. In the last section, our numerical results are given, 
followed by a discussion. 



II. FORMALISM FOR MESON-BARYON INTERACTION 
A. Lagrangian and Feynman diagrams 

We consider the PB — > PB and VB — > VB interaction by exchanging a vector meson. 
The corresponding Feynman diagrams are shown in the FigfTJ 




(a) (b) 



FIG. 1: Feynman diagrams for the pseudoscalar-baryon (a) or vector- baryon (b) interaction via 
the exchange of a vector meson (Pi, P2 are D~ , D° or D~ , and V\, V2 are D*~, D*° or D*~, and 
B\, B2 are S c , A+, H c , Ef c or Sl c , and V* is p, K*, (f> or oS). 



In order to evaluate these Feynman diagrams, we give the three types of vertices for BBV, 
PPV and VVV interaction from 9|. The Lagrangians for the interaction of vector mesons 
between themselves (three - vector vertex), pseudoscalar mesons with vectors and baryons 
with vectors are: 

Cvvv = tg(V»[V\dM) 
Cppv = -ig(V»[P,d,P}) 

Cbbv = g{{Bj,[V»,B}) + (B lfl B)(V»)) (1) 



where B and P are the standard matrices including the pseudoscalar and baryon nonets 
27|. g — My/2f is the coupling used in the hidden gauge with the pion decay constant 
/ = 93 MeV and the mass of the light vector meson taken as My = 770 MeV. The g fulfills 



the KSFR rule 
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which is tied to vector meson dominance 



When we go to SU(4) 



we can still use the Lagrangian for VPP of Eq. (CO) and the V and P matrices extended to 
SU(4): 



2I_ _i_ J?a_ _i_ Vc 
V2 ~ l ~ V6 ^ VT2 



and 



7T 



K~ 



7T 



V2 ^ ^ 



I Vc 

12 ^ V4 



K c 



if 



-2r? 8 



+ 



+ 



D c 



D~ 



D: 



_3ijc_ _|_ g. 



(2) 



£ |_ U8. _|_ 



v 7 ^ V6 Vl2 \fl 



_|_ j£c. 



V 



D 



*o 



V2 -v/6 



*o 



*0 



A' 



D*~ 



*0 



D 



D* 



D*~ 



3ui c 
\f\2 



• (3) 



Let us recall that here fj c stands for the SU(3) singlet of the 15th SU(4) representation 
and we denote fj' c for the singlet of SU(4). On the other hand, u!$ plays the role of the r] 8 
for the vectors, while u) c the role of fi C) and we denote by Cj' c the SU(4) singlet. We take n°, 
i], rj', rjc as a basis for the neutral pseudoscalar mesons, where r( is the singlet in SU(3), 
(uu + dd + ss) I \/3, and r) c stand for cc. Recalling the standard quark composition of the 
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SU(4) mesons 







n° = 


V ^ 






V8 = 


71° 




i 


Vc = 


y/12 






€ = 


71° 



uu — dd) 
uu + dd — 2ss) 
uu + dd + ss — 3cc) 
uu + dd + ss + cc) , (4) 



we find 



Vs = V 

11 = 2 Vc + ~Y Vc 



h-Vz^ + t), (5) 



Vc 2 

in the physical basis. On the other hand, for vectors we use the physical basis p°, u>, <fi and 
J/ip, where 

p° = —= (uu — dd) 
v2 

uj = —j= (uu + dd) 

= ss 

J/ip = cc , (6) 



which can be written in terms of u 8 , u c and u' c as 1 



u = l(V6uj c + 2V3uj 8 + 3V2uj' c ) 
o 

= -(V3u c - 2v / 6a; 8 + 3u' c ) 

= l(-Vau e + u' e ) . (7) 

The use of Lagrangians to give the BBV vertex in SU(4) is more cumbersome than in SU(3) 
and thus it is simpler to use SU(4) Clebsch Gordan coefficients. Yet, this requires a certain 
phase convention for the physical states with respect to the isospin states implicit in the 



[1] Latter on, in order to use the SU(4) Clebsch Gordan coefficients we shall change a phase to the fj c and 
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SU(4) tables, which makes convenient to use the same procedure to evaluate the PPV 
vertex. 

In the PPV vertex we go from the 15 <S) 15 representation of pseudoscalars to the 15 
representation of vectors. Yet, the nature of the couplings (with the explicit commutator) 
has as a consequence that only the 15 p (antisymmetric) representation for the vectors is 
needed (one can alternatively use explicitely the 15f and the 15d representations and the 
15d contribution vanishes at the end 2 ). The resulting t amplitude for P1P2 — > V is given by 

tp 1 p 2 v = 9i5 F C 15f (15 <g> 15) (q 1 + q 2 ) ■ e , (8) 

where q\ and q 2 are the four-momentum of the initial and final pseudoscalar mesons respec- 
tively, and Ci5 F (15(g) 15) is the SU(4) Clebsch Gordan Coefficient that we take from [28| and 
gis F is the reduced matrix element that by comparison with the result of the Lagrangian is 
given by 

g 15p = -2V2g . (9) 

However, the use of the SU(4) tables requires a phase convention. We find a compatible 
and convenient phase convention of the isospin states implicit in the SU(4) tables and those 
used by us in Eqs. © and (jHJ) by means of: 

|tf°} = -|l/2,-l/2> ; |tt+) = -|1,1) ; = -\1,0) ; 

\D+) = -|0,0> ; \D°) = -11/2,1/2) ; \fj c ) = -|0,0) ; 

and equivalent phases for the corresponding vectors, K*°, p + , p°, D* + , D*° and u c . The 
necessity for the change in phases stems from demanding that the 15 <E> 15 — > 1 combination 
of SU(4) isospin states is a symmetrical expression in the physical states {2^. The use of 
this convention (and also the convention for baryons that we give later) leads to the same 
amplitudes in charge basis given by the Lagrangians of Eq. ([1]) with the P and B matrices 
written in the SU(3) basis. 

When we go to the BBV vertex (we look for BB — > V), we need now the three repre- 
sentations, 15i, 152 and 1, and we must note that when the 8 representation of SU(3) is 
involved, only the F coefficients are needed. In this case we have 20' <g> 20 ->■ 15i, 15 2 , 1, and 

[2] We thank J. Nieves for pointing this to us. 
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the t amplitude for the BBV vertex is given by 

t Bl S 2 v = {915, Ci 5l (20 / ® 20') + g 152 C 152 ( 20' ® 20') + <?i C 1 (20 / <g 20')}u r/ {p 2 ) 1 ■ eu r { Pl ) . 

(10) 

Once again by writing the expression for 20' <g> 20 -> 1 in terms of the SU(4) isospin states, 
and demanding that the expression is symmetrical in the physical baryons, we obtain a 
convention of phases. The one we have chosen, partly motivated to agree formally with 
earlier SU(3) results, is given by changing the phases of the states 

|S~> = -(1/2, -1/2) , |fl+> = -|0, 0) , |Sg> = -|l/2, -1/2) , |5») = —11/2, -1/2) , 

|A c -) = -|0,0) , |£+) = -|l,0) , |E++) = -|l,l) , |E--) = -|l,-l) , 
\n) = -|l/2, -1/2) , |H°) = —11/2, -1/2) , |E+) = -|1, 1) , |E+) = -|1, 1) , 

|E°) = -|1,0) , |S°> 11,0> . 

The reduced matrix elements of Eq. ffTOl . g l5l1 g l52 and g± are evaluated demanding: 

1) The coupling pp — >■ J/ip should be zero by OZI rules, 

2) The coupling pp — > <p should be zero by OZI rules, 

3) The coupling pp — > p° should be the one obtained in SU(3). 
We finally obtain 

015! =-01 9i5 2 = ^V3g; gi = 3V5 g . (11) 

with g = My/2/ and / = 93MeV" the pion decay constant. 

The diagram of Fig. [1] (a) requires the exchange of the vector meson with the two vertices 
given by Eqs. (jHJ) and (TlOl) . In the sum of polarizations in the vector meson exchanged, 

X} e M e " = -0A** + Tsjr ■ ( 12 ) 



A 



V 



We can keep just the /i = v — component since we assume that the three momenta of the 
particles are small compared to their masses. Similarly, the q 2 /My term in the vector meson 



propagator is neglected (further on, when we consider the transitions from heavy mesons to 
light ones, we perform the exact calculation). The transition potential corresponding to the 
diagrams of Fig. [Tj are given by 

Vab(P 1 B 1 ^P 2 B 2 ) = jjf0?l+<?2); ( 13 ) 
Vab(ViB 1 ^V 2 B 2 ) = Tf|(9i + 92)?L ' ^2 • (14) 



Where the indices a, b stand for different groups of Pi(Vi)Bi and ^2(^2)^2, respectively. 
The ql, q® are the energies of the initial, final meson. We list the value of the C a b coefficients 
for different states of isospin, I, and strangeness, S in the Appendix. Here we study six 
different cases with (I, S) = (3/2, 0), (1/2, 0), (1/2, -2), (1, -1), (0, -1), (0, -3). 

B. The G function and the unitary T amplitudes 

The G function is a loop function of a meson (P) and a baryon(5) which we calculate in 
dimensional regularization by means of the formula 



<Pq 1 i 

G (m = z2M B I {27TriP _ q) 2 1 _ M 2 +l£ g2 _ M 2 + - (15) 



2M B r , Ml Ml-Ml + s, Ml 
*-{a At + ln^+ p B In ' 



16tt 2 1 » n 2 2s M 2 B 

s ■ 

ln(-s - (Ml - M 2 P ) + 2q^~s) - ln(-s + (M 2 B - M 2 P ) + 2q^~s)] } , (16) 



+ [ln(s - (Ml - Ml) + 2q^) + \n(s + (M 2 B - M 2 P ) + 2q^fs) 
V s 



where 



s = P\ (17) 

_ = y/ (.- {UB+Ur m.-( UB - UpY) with imfe) > o (i8) 

2y s 

In Eq. ffTBI) . q is the four-momentum of the meson, and P is the total four-momentum of 
the meson and the baryon. The \x is a regularization scale, which we put 1000 MeV, and 
is of the order of —2, which is the natural value of the subtraction constant 3j. When 
we look for poles in the second Riemann sheet, we must change q by — q when -y/s is above 
the threshold in Eq. (jTB"]) [30I . See further comments regarding the subtraction constant in 
Subsection D. 



S 



Here we also regularize the G loop function in a different way by putting a cutoff in the 
three-momentum. The formula is: 

d 4 q 1 1 



G (PiB) = i2M B 



(2tt) 4 (P - q) 2 -M 2 +ieq 2 -M 2 + is 
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g 2 <ig 2M B (wp+wg) 
4vr 2 ujpu B ((P ) 2 - (u P + u B ) 2 + ie) 



(19) 



where 



w P = Jq 2 + M 2 



P 5 



WB = 0f 2 + M 2 , (20) 

and A is the cutoff parameter in the three-momentum of the function loop. 

For these two types of G function, the free parameters are a M in Eq. ( TT6|) and A in Eq. 
f fT9|) . When we choose a M or A, the shapes of these two functions are almost the same close 
to threshold and they take the same value at threshold. 

Then we can get the unitary T amplitudes by solving the coupled channels Bethe-Salpeter 
equation in the on shell factorization approach of js], [mI, 32] 

T =[1-VG}~ 1 V . (21) 

When we look for poles in the complex plane of y/s, poles in the T matrix that appear 
in the first Riemann sheet below threshold are considered as bound states whereas those 
located in the second Riemann sheet and above the threshold of some channel are identified 
as resonances. 



C. The coupling constant and the width of the poles 

From the T matrix we can find the pole positions zr. In this work, we find all of these 
poles in the real axes below threshold, in a few words, they are bound states. In view of 
that, for these cases the coupling constants are obtained from the amplitudes in the real 
axis. These amplitudes behave close to the pole as: 

T ab = . (22) 

We can get the coupling constant as: 

g 2 a = lim (T aa x (V7s-z R )). (23) 



This expression allows us to determine the value of g a , except by a global phase. Then, the 
other couplings are derived from 

g b = hm (^*) . (24) 
As all the states that we find have zero width, we should take into account some decay 



mechanisms. Thus, we consider the decay o 



means of box diagrams as it was done in 33 



the states to light meson - light baryon by 



34j . The Feynman diagrams for these decays 



x- 



are shown in Fig. [2J We assume that P3, V3 and B 3 are on-shell and neglect the three - 
momentum of the initial and final particles. Then, using Eq. (pQ), the transition potential of 
these diagrams can be written as: 

T/ C ac C cb M v , n (y^ + M B3 ) 2 -M|, 3 

Vacb(P 1 B 1 -+P 3 B 3 -+P 2 B 2 ) ~ YqJI X Wft.Bs) X A^/sMb 

2E Pl + {M B , - M Bl ){M Pi + M v , - M%)/M*. 
M\ + M 2 Ps - 2E P3 E Pl - My* 
-2E P2 + {M B3 - M B2 ){M\ + Ml, - M|) /My, 

M P2 +M P3 -2E Ps E P2 -My, ' { } 

and the same for vectors (see Fig. [2j (b)) changing E Pl , E P2 , E P . A by Ey x , Ey 2 , Ey z and M Pl , 
M P2 , M Pi by M Vl , My 2 , My 3 , respectively. Here c stands for a different group of P 3 (V 3 )B 3 . 
Then, the kernel V in the Bethe Salpeter equation, Eq. fT2TT) . becomes now: 

V a b(P 1 B 1 ^P 2 B 2 ) = ^p( E Pi + Eft) + X] Vacb ' ( 26 ) 

and similarly for the VB system. In Eq. (|25|) we have factorized the two P\B\ — > P3-B3 
and -P3-B3 — > P2B2 transition amplitudes outside the loop integral by taking their values 
when the system P3-B3 is set on-shell. This is a good approximation, exact for the imaginary 
part of the diagram, which is our main concern, since we are interested in the contribution 
of these diagrams to the width of the resonances. The loop integral only affects then the 
P 3 , P 3 propagators leading to the same G function defined in Eq. ([16]) . Since the on-shell 
mass of the intermediate states is far away from the energies investigated, ReG(P3, P3) is 
small and we have checked that it is sufficiently smaller than the tree level contribution 
from the diagrams of Fig. ([T]), such that it can be ignored. For example, V^ p , t , c ^ttN-^dy, c ) = 
(0.38 + 2.9i) GeV- 1 at the N* pole position with yfe = 4.265 GeV. 

Further on, we will include the r] c N, t] c A channels for PB — > PB, and J/ipN, J/ipA for 
VB — y VB in the calculation. 
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FIG. 2: The Feynman diagrams of pseudoscalar-baryon (a) or vector-baryon (b) interaction via a 
box diagram. Pi, P2, V±, V 2 , B\, B2 are the same particles than in Fig. [TJ P3, V3 and B3 are light 
particles belonging to the SU(3) octet of pseudoscalar mesons, vector mesons and stable baryons, 
respectively, and V*, V 2 * are D* or D*. 

D. Discussion about the use of the SU(4) flavor symmetric Lagrangian 

Once the formalism has been exposed we would like to make some comments to justify 
the approach. While the hidden local gauge approach is well settled with SU(3) flavor, its 
extension to SU(4) is not quite justified. The local hidden gauge theory in QCD is based on 
spontaneous symmetry breaking of chiral symmetry, which is not expected to hold in the 
charmed sector, and even if this were the case, the breaking pattern would be masked by the 
large charm mass. In view of this, the approach followed requires some justification. The 
first thing we must bare in mind is that the large mass of the charmed quark is to be blamed 
for the lack of symmetry. Hence, if hadron dynamics still has some traces of SU(4) symmetry 
it should be in particular vertices or amplitudes not tied to the quark mass. In this respect 
it is interesting to mention that SU(4) symmetry works fairly well for the vertices VP7 (or 
equivalently the VVP) involved in the radiative decay of vector mesons with charm to a 
pseudoscalar and a photon [35|]. The agreement is as good as in SU(3). Given the analogies 
between the VVP and VVV vertices provided by heavy quark symmetry, it is fair to think 
that using SU(4) symmetry to evaluate the VVV and VPP vertices would be a fair starting 
point. Similarly we could also assume the same symmetry to hold in the VBB vertices. 

One should also note that in the case of meson meson interaction the present approach 
provides the same results, up to a mass term of no practical consequences, as the heavy quark 
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formalism used in |36H38| for the case of interaction of light mesons with heavy meson. From 
the perspective of our approach, we could rephrase it by stating that these two approaches 
provide the same VVV (or VPP) vertices with two heavy vectors and a light one. For 
the evaluation of the width in the present approach one is using these vertices, with the 



only difference with respect to 36M38j that one of the heavy vectors is exchanged in the 



t— channel in our approach, while the two heavy vectors (or pseudoscalars) were external 



particles in 



3614381] . Only vertices involving three heavy vectors (D*D*J/ip) would require 



the extra help of SU(4) symmetry. Thus, most of the information used is supported by 
phenomenology and other approaches. 

Furthermore, we should bare in mind that the largest fraction of the results that we ob- 
tain, concerning couplings to different channels, is based on SU(3), since we can relate these 
channels through SU(3). Then, implementing SU(4) symmetry as we do, automatically ac- 
complishes this. Only when we give a jump to another charm sector the SU(4) symmetry 
would play a role and there we would invoke the arguments used above to support it. Cer- 
tainly, when these vertices are used in Feynmann diagrams and the masses of the exchanged 
vectors are very different, the approximate SU(4) symmetry that we had in the vertices will 
be badly broken in the amplitudes. This happens also in SU(3): The vertices are manifestly 
SU(3) invariant, but when SU(3) is broken in the amplitudes because of the different masses 
of the particles belonging to the same multiplet (for instance in the unitarization procedure), 
the underlying SU(3) symmetry is broken and two octets that were degenerate in the exact 
SU(3) limit give rise to two different states in the strangeness S = — 1 sector: one of the two 
A(1405) states and the A(1670). 

Yet, one should be ready to accept larger uncertainties than in SU(3) and allow some 
fitting freedom in the approach. This can be done by means of the subtraction constants 
of the G function, that effectively tune the strength of the potentials that one is using in 
the approach. This also means that the natural values of these constants should only be 
used as indicative and then a real fit to the data should be done, which cannot be done in 
the present case since we have no experimental data. However, one can rely on previous 
work along these lines in which several groups have done this work and provide the new 



scale o: 

of m 



re subtraction constants to be used in the charm sector. In this sense, the works 



39M42l| . choosing these parameters to reproduce properties of known resonances like 



the D° (2317), the X(3872), the L>*(2460) and £>*(2640), and A c (2593), have given us the 
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scale for these subtractions constants that we use here. 



III. RESULT AND DISCUSSION 



A. The pole positions and coupling constants 

Here we show the results for the different sectors. By using the two G functions of Eqs. 
ffTB]) and ( TO?]) , the poles appear in both cases below threshold in the first Riemann sheet and 
therefore they are bound states. We show the pole positions for different values of a (A) in 
Tables [I] and [Til 

We take a range of values of a, or accordingly the cut off, in line with values used in 
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39h41| and we find six poles in our calculation. The uncertainties in the pole positions 
in the case of the first and third poles for both PB and VB systems, are of the order of 100 
MeV, which are typical in any hadron model. These two poles are rather stable. However, 
for the second state, the uncertainties are much larger and the pole position is very unstable. 



For the discussions that follow we choose an intermediate value of a, which we take 
a = —2.3, to study the nature of these poles in detail. In Tables UTTI and ITVj the values of 
the coupling constants are listed by using Eqs. ( 1231) and ( 124T) . From Table UTTI we see that 
both the A* (4269) and the A* (4403) depend on one channel, DE C and D"E' C , respectively. 
These two states are both stable as we can see in Table [TJ In contrast, the A* (4213) depend 
on two channels, D S A^ and DS C . The mass of this state changes appreciably by using 
different values of the free parameters (a or A). 



B. The decay widths of these states to light meson - light baryon channels 

These states decay to two different types of channels, one is the light meson - light 
baryon channel, while the other is the cc meson - baryon channel. For the VB states, there 
is another possibility to decay into PB channels, for instance, D*B — > DB. The analogous 
decay channels in the VV — > VV hidden charm sector driven by pseudoscalar exchange were 
studied in 4^] and found to be extremely small because of the small phase space available. 
Analogously, the terms involving a vector exchange contains an anomalous VVP vertex and 
were also found very small in 33j. Hence, we do not consider them here. In this subsection 
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(I,S) 


a = -2.2(A = 0.7 GeV) 


a = -2.3(A = 0.8 GeV) 


a = -2.4(A = 0.9 GeV) 




ZR 


ZR 


ZR 


(1/2,0) 


4291(4273) 


4269(4236) 


4240(4187) 


(0,-1) 


4247(4120) 


4213(4023) 


4170(3903) 




4422(4394) 


4403(4357) 


4376(4308) 



TABLE I: Pole position from PB — > PB using the two different G functions of Eqs. (fTUj) and (fTUj) , 



The units are in MeV. 



(I,S) 


a = -2.2(A = 0.7 GeV) 


a = -2.3(A = 0.8 GeV) 


a = -2.4(A = 0.9 GeV) 




ZR 


ZR 


ZR 


(1/2,0) 


4438(4410) 


4418(4372) 


4391(4320) 


(o,-i) 


4399(4256) 


4370(4155) 


4330(4030) 




4568(4532) 


4550(4493) 


4526(4441) 



TABLE II: Pole position from VB —> VB using the two different G functions of Eqs. (116)) and 
([15]) . The units are in MeV. 



we only consider the decay of these states to the light meson - light baryon channel as 
depicted in the Feynman diagrams of Fig. |2j These diagrams provide a negligible real part 
compared to the tree level potentials. The imaginary part gives rise to a width of the states. 
Hence, we only consider the effect of this box diagram on the states found before. 

In Figs. EJandlU we show the results of \Ta\ 2 as a function of yfs for the different channels, 
and we list their decay widths to the different channels for all the sectors in Tables [V] and 
IVII From these pictures and tables, we find that the six states are all above 4200 MeV. 
However, their widths are quite small. In principle, one might think that the width of these 
massive objects should be large because there are many channels open and there is much 
phase space for decay. However, it is difficult for the cc components to decay to the uu, dd 
and ss ones, something that within our model is tied to the necessity of the exchange of a 
heavy - vector meson. Note that the pole positions are obtained without including the box 
diagrams by extrapolating to the complex plane. The inclusion of the box diagram renders 
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(I, S) 


zr (MeV) 




9a 




(1/2,0) 




DY, C 


DA+ 






4269 


2.85 







(0,-1) 




D S A+ 


DE C 


DE' C 




4213 


1.37 


3.25 







4403 








2.64 



TABLE III: Pole positions, zr and coupling constants, g a , for the states from PB — > PB. 





zr (MeV) 




9a 




(1/2,0) 






D*h + 






4418 


2.75 







(o,-i) 






D*E C 


D*E' 




4370 


1.23 


3.14 







4550 








2.53 



TABLE IV: Pole position and coupling constants for the bound states from VB — > VB. 



this extrapolation more difficult, and thus we obtain the width of the states by plotting 
\T\ 2 versus the energy with T obtained in the real axis including the box diagrams. The 
individual partial decay widths are obtained including one by one the different box diagrams. 



C. Decay width to cc meson - light baryon channels 

In this subsection we discuss the decay width of these states to cc meson and light - baryon 
channels. The three states from the VB system decay to J/ipN . The decay of these VB 
states to rj c N is also possible by means of a BBP vertex (exchange of a pseudoscalar meson) 
but as we will see in the Subsection IV. B this vertex is very small. We could also have 
this decay exchanging a vector meson instead of a pseudoscalar one, but then the amplitude 



would contain an anomalous VVP vertex, which is also very small [33j. Similarly, the decay 
width of the PB states to the VB channels must be very small because of the same reasons. 
We will consider their decay to J/ipN in Section IV. B and we anticipate that this decay 
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FIG. 3: \Tu\ 2 for the different channels in the (I = 1/2, 5 = 0) sector including the box diagrams. 
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FIG. 4: \Tu\ 2 for the different channels in the (I = 0,S = —1) sector including the box diagrams. 
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M 


r 












(1/2,0) 










r)N 


7]'N 








4269 


4267 


34.3 


3.8 


8.1 


3.9 


17.0 




(0,-1) 








KN 




r)A 


rj'A 


if- 




4213 


4213 


26.4 


15.8 


2.9 


3.2 


1.7 


2.4 




4403 


4402 


28.2 





10.6 


7.1 


3.3 


5.8 



TABLE V: Pole position (zr), mass (M), total width (r), and the decay width for each particular 
light meson - light baryon channel (Tj) for the states from PB — > PB. The units are in MeV. 



(I, S) zr Real axis Tj 

m r 



(1/2,0) 


4418 


4416 


28.4 


pN 
3.2 


loN 
10.4 


K*T, 
13.7 






(o,-i) 








K*N 


P S 


loA 


cj)A 


K*E 




4370 


4371 


23.3 


13.9 


3.1 


0.3 


4.0 


1.8 




4550 


4549 


23.7 





8.8 


9.1 





5.0 



TABLE VI: Pole position (zr), mass (M), total width (r), and the decay width for each particular 
light meson - light baryon channel (Tj) for the states from PB — > PB. The units are in MeV. 

width is very small. For these reasons, we only consider the J/tpN, J/ipA channels for the 
V B states, and i] c N, i] c A channels in the case of states from the PB system. Thus, these 
new channels are added to the previous calculation in the Subsections III. A and B. 

The pole positions of these states only change a bit compared to those given in the 
Subsection III. A, since the potentials from these channels are much smaller. Nevertheless, 
these channels provide some extra width because, in spite of the smaller phase space for 
the decay, the three momentum transfer in the propagator of the D*(D*) exchange is much 
smaller than in the case of transition to light meson - light baryon channels. The transition 
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potential becomes: 



C b9 2 

Vab(PB-+ Vc B) = - 2 [_ M2 (E P + E Vc ) , (27) 
Pd* m D* 



where 

p 2 D , = Ml +M 2 P - 2E Vc E P , (28) 

and similarly for the VB system but changing p D *, M D *, E P and E Vc by p D , M D , E v and 
£ju respectively. Here we also neglect the three-momentum of the final and initial particles 
because we consider energies close to the threshold. We list the results in Tables I VI It IVIII1 
[IX] and |Xl We observe that the coupling constants change a bit, but what is more relevant 
is that these new channels give an extra contribution to the width, smaller, but of the same 
order as the one obtained previously. The relatively large decay width to the rj c N channel is 
a good feature with respect to the possible observation of these resonances since there will 
be less background in rj c N than in ttN, t]N, KT,, the observation of the resonance in the 
r] c N channel could be favoured. 

In Tables IVHI and IVIII1 the pole positions are obtained without the box diagrams, but 
including the 1] C N, rj c A channels. Now the pole positions becomes complex because the 
new channels are open. We can see that the partial decay width into these channels is 
approximately twice the imaginary part of the pole position. The total widths are again 
obtained by looking at the width of |T| 2 in the real axis when the box diagrams are includec . 



18j |. which has been corrected in |19l. |20| . 



We would like to mention that in the approach of 
some hidden charm states are also found, bound by about 1000 MeV. It is not easy to 
understand such a large binding on physical grounds, which is not supported in any case by 
the strength of the potentials. 



IV. PRODUCTION CROSS SECTION IN pp COLLISIONS 

A. Estimate of the pp — > N*^ (4265)p cross section 

We shall estimate the production cross section of these resonances at FAIR. With a p 
beam of 15 GeV/c one has y/s = 5470 MeV, which allows one to observe resonances in 
pX production up to a mass Mx — 4538 MeV. We shall make some rough estimate of the 
cross section for the pp — > pN*^ production for the C = 0,5 = resonances that we have 
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(7,5) z R (MeV) 



(1/2,0) 7J£ C DA+ r] c N 



4265 - 11.6* 


2.96 - 0.21?' 


-0.08 + 0.06? 


-0.94 + 0.03? 






2 Q7 

— ■ -J \ 




D Q4 




(o,-i) 


D S A+ 




DE> C 


%A 


4210 - 2.9? 


1.42 - 0.03? 


3.28 - 0.002? 


-0.15 + 0.13? 


0.57 + 0.04? 




1.42 


3.28 


0.19 


0.57 


4398 - 8.0? 


0.01 + 0.004? 


0.06- 0.02? 


2.75 - 0.15t 


-0.73 - 0.07i 




0.01 


0.06 


2.75 


0.74 


TABLE VII: Pole position, zr and coupling constants, g a , to various channels for the states from 


PB —7- PB including the r/ c N and ry c A channel. 








(7,5) 


z R (MeV) 


Real axis 


Ti 








M r 






(1/2,0) 






r] c N 






4265 - 11.6? 


4261 56.9 


23.4 




(o,-i) 






r] c A 






4210 - 2.9? 


4209 32.4 


5.8 






4398 - 8.0? 


4394 43.3 


16.3 





TABLE VIII: Pole position (zr), mass (M), total width (r, including the contribution from the 
light meson and baryon channel) and the decay widths for the 77 C A^ and r/ c A channels (Ti). The 
unit are in MeV 



obtained from the pseudoscalar baryon interaction. Since one important decay channel of 
the N* 5 is ttN, we evaluate the cross section for the mechanism depicted in the Feynman 
diagram of Fig. [5j 

The coupling of the N* 5 — > ir°p is obtained projecting over 7r p the isospin state 7 = 1/2, 
which provides the isospin coefficient Ci = a/1/3. The coupling N* 5 — > ttN we get from the 
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(/, S) z R 



(1/2,0) 


D*T, C 


D*A+ 


J/^pN 




4415 - 9.5* 


2.83-0.19? 


-0.07 + 0.05*' 


-0.85 + 0.02*' 






2.83 


0.08 


0.85 




(0,-1) 

V"> *-) 


D*A + 


D*E r 

1 1 — 'C 


D*E' 


J/ipA 


4368 - 2.8* 


1.27-0.04* 


3.16 - 0.02*' 


-0.10 + 0.13* 


0.47 + 0.04*' 




1.27 


3.16 


0.16 


0.47 


4547 - 6.4* 


0.01 + 0.004* 


0.05 - 0.02*' 


2.61 - 0.13* 


-0.61 - 0.06*' 




0.01 


0.05 


2.61 


0.61 


TABLE IX: Pole position (z R ) 


• and coupling constants (g a ) to various channels for the states from 


PB — > PB including the J/tpN and J/ipA channels. 






(I,S) 


ZR 


Real axis 


Ti 








M r 






(1/2,0) 






J/ipN 






4415 - 9.5* 


4412 47.3 


19.2 




(o,-i) 






J/ipA 






4368 - 2.8* 


4368 28.0 


5.4 






4547 - 6.4*' 


4544 36.6 


13.8 





TABLE X: Pole position (zr), mass (M), total width (r, including the contribution from the light 
meson and baryon channel) and the decay widths for the J/ipN and J/ipA channels (Tj). The unit 
are in MeV 



partial decay width of the N* 5 into this channel, T nN 

2 27rMjv*_r 7r jv , , 

with p° n = A 1 / 2 (M^ r « _, m 2 , M^)/2M N *_, the value of the on-shell pion momentum from the 
N* 5 — > irN decay. By taking the standard irNN vertex, V wNN = ig^r x (g n ~ 13), we 
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p — « — 1 « p 

FIG. 5: The pp — > N*^~p mechanism. 



obtain 

dv^N-jp = glM\ Y^C} 2p.p'-2M 2 p' 

dcos6 4 s p° n (2M 2 - y/sE(p') + 2p.p') 2 p 1 ' 

where p, p' are the initial, final momenta of the p in the center of mass frame ( of the order of 

2570, 620 MeV/c for Mx — 4300 MeV). The biggest cross section corresponds to the forward 

p direction, which is the most indicated for the search. If we are interested in searching for 

these resonances, looking for p forward is the most recommendable measurement and one 

should look for a bump into the da/dcosOdM 2 magnitude, where Mj is the invariant mass of 

the 7rN coming from the decay of the produced iV* g + state. Assuming a Lorentzian shape for 

this resonance, with total width T w *+, we would obtain at the peak of the nN distribution 

dcosOdMf ~ 7rM w .+r tot dcosO 

cc 

which leads to the following cross section: 0.13 fib/GeV 2 for iV* g + (4265). 

In the above calculation, we did not consider the form factor for the nNN vertex. The 
form factor is: 

A 2 — rn 2 

Fpp7T = \2 _ n 2 ■ ( 32 ) 
with the A n = 1.3GeV. We can multiply by F 2 p7T the cross section in the Eq. (|3T|) and we 
find about 0.05 fxb/GeV 2 . 

Because in such high energy transfer reaction the one-pion exchange with the monopole 
off-shell form factor of Eq. ( 1321) may not be a good approximation, here we also make a 
calculation with the Reggeon exchange. Using a Reggeon propagator R n (s,t) \^ instead 
of the usual pion propagator. Then the Eq. fl30|) becomes 



= f ~^T' E ^l - ^«** - 2Ml)\Ms,t)f-, (33) 
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where 



^ ^ sin(|a^(t))r(-^ + 1) 

a w (t) = -0.015 + 0.72*, (35) 

a' n (t) = 0.72, (36) 

t = 2M 2 N - yfsE(p') + 2pp'cos6, (37) 



with the slope parameter 



R w ~ 2 t . From Ref. 



in the units of (GeU 2 ). When i — > m 2 , we can get = and 



44j, the order of So is about 2 — 20GeV~ 2 . To narrow down its range, 
we use the information of pp collision with y/s < 3 GeV where the one-pion exchange can 
reproduce experimental data reasonably well [451 ] . Demanding the Reggeon propagator to 
give similar results as the usual ir propagator for y/s < 3 GeV, we have s — 5 — 10GeV~ 2 . 
Then by using Reggeon propagator for the maximum PANDA energy y/s < 5.47 GeV, we 
get the cross section to be about 0.006 ~ 0.017 pb/GeV 2 corresponding to s = 5 ~ 10 
GeV -2 . This is about a factor 3 ~ 9 smaller than the result by one-pion exchange. 
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FIG. 6: The different Feynman diagrams of the reaction pp — > ppi] c 



Then we can estimate the cross section of pp — > ppr] c . The different Feynman diagrams 
for this reaction are shown in the Fig. [6J Using Eq. (131 j) and T^ cP of the resonance instead of 
r n N we can obtain the differential cross section at the peak of the resonance, corresponding 
to the resonant mechanism of Fig. [6]a), and it is about 0.8 pb/GeV 2 without form factor and 
0.3 [ib I 'GeV 2 with the form factor, and 0.04 ~ 0.10 fib/GeV 2 with the Reggeon propagator. 
This magnitude is of about the same order of magnitude as typical cross sections measured 
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for da/dcos8dMf in the pd — > pdn ?: or pp — > d7c + 7c° reaction 46|, |47|. In order to see the 
role played by the hidden charm resonance in this process we can compare it with the cross 
section coming from a standard mechanism of Fig. E](c,d). The vertex of ppr\ c is used by 

£vcpp = 9vcPpU P l^ 5 d^ Vc Vp, (38) 
where g VcP p can be calculated from the reaction r\ c — > pp by 



9w - v ~ww ( ' 

where the p° n = A 1 / 2 (m^ c , M^, M?)/2m r?c the value of the on-shell p momentum from the 
r] c — > pp decay. And the width T Vc = 26.7MeV and the branch ratio Br^pp = 1.3 x 10 -3 
are both from PDG. The form factor of the vertex NNtt is also used Eq. fj32|) . We also add 
the form factors for N* 5 and p exchange in the Figj6j 

Here A p = A N = 0.8GeV. 




FIG. 7: The total cross section vs the beam momentum of p for pp — > ppr] c . The solid and dashed 
lines are calculated by one-pion exchange without and with form factors, respectively. The dot- 
dashed and dotted lines are results by Reggeon propagator with sq = 5 and 10 GeV -2 , respectively. 
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Through the calculation, the contributions from Fig. [6] (c), (d) are very small, almost 
10~ 4 /j,b, the main contribution comes from the N* 5 . The total cross section is about 0.0029/^6, 
0.013/ife, 0.072/Z& and 0.71/ife for a p beam of 15 GeV/c as shown in FigJ3, corresponding to 
the Reggeon propagator with sq = 5 GeV -2 and Sq = 10 GeV -2 , the usual tt propagator 
with and without form factors. Note that the integrated cross section involves finite angles, 
rather than zero in the forward direction considered before, where the effect of the form 
factor is more important. The Dalitz plot, the invariant mass spectrum of prj c , pr] c and pp 
are all shown in Fig. [8] where the peaks of N* (4269) are very clear. 




15.51616.517 17.51818.519 15.51616.51717.51818.519 



(c) Mp^GeV 2 ) (d) M^ c (GeV 2 ) 

FIG. 8: The Dalitz plot(a), the invariant mass spectrum of pp(b), PVc( c ) &nd pr] c (d) for the reaction 
pp —7- ppr] c at the beam momentum of p being 14.00GeV at lab system. 
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B. J/ip production in pp — > ppj/tp. 



Another estimate that we want to do is the cross section for J/ip production in the 
pp — > ppJ/ip reaction around the energy of the iV*(4265) excitation. We use again Eq. f |5Tj) 
but we need to evaluate Tjm p . This requires a different formalism to the one used so far. 



The mechanism for R — > J/ipp is obtained by analogy to the work done in 33j, |34j where 
the transition from vector - vector to pseudoscalar - pseudoscalar states is done. Concretely, 
given the fact that the A^"(4265) is basically a -DS C molecule in our approach, we obtain 
the coupling of the resonance iV* g + (4265) to J/ipp through the diagram of Fig. |9j 



£ + „(P-q) 




DO (q y--KAAAAWj/ ¥ (p J/¥ ) 



FIG. 9: pJ/tp going to the resonance iV* g + (4265). 



This diagram requires the coupling of iV* 5 + (2465) to the -DS C state in I — 1/2, and the 
transition J/ipp — > DX C which is mediated by the D meson that comes from the coupling of 



J ftp to DD. The diagram also involves the DNY* C coupling which has been studied in 48]. 
The J/ip —t- DD coupling can be obtained from the Lagrangian 

Cp PV = -ig(V< i [P,d tl P)) , (42) 
used in Section II, with g = My/2/ and / = 93 MeV, which leads to 

- ttj/^DD = i 2g ■ (43) 
The vertex DNH C is obtained from 48j and has the form 



„/0 jj p 

- iV D o pSt = 3 ■ ^(1 - ^^J~ ( 44 ) 
with (3 = 1 and q'°, q', the incoming energy, momentum of the D meson and M' the mass 



49-511. Thus, 



of the S c . For D and F we take the standard values D = 0.8 and F = 0.46 

- rf DO P s+ = ■ Q (45) 
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We need the 7 = 1/2 state of -DS C given with our phase convention by 

|D£ C ; 1/2, 1/2) = + • ( 46 ) 

The other possible vertex, the _D + p£+ + vertex, is \[2 times the D°pY<^ one. With all these 
ingredients one obtains 

r- f d*q 0.26^ _ _ M Sc 1 



(2tt) 4 2f * H E^ c {q) q 2 -rn? D + ie 
X (pj - qf -ml + t ePO-qO- fyTfa) + ' (4?) 

where we use a form factor F(q) = with A = 1.05 GeV [48] in the integral of Eq. (T47j). 

Upon neglecting the small three momenta p.j/^ compared to the J/ip mass and performing 
the q° integral, Eq. (|47|) can be written as 

1 0.26 _ _ f d 3 q ^ 2 M Sc 1 1 1 



-itj/^ R = -=— — ga-e / j^—^q 



^3 f J (2n)^ E Se {q)2u D (q)p j + 2u D (q)p°j-2u D (q) 

1 1 



x 



P o _ p o _ uM _ E ^ q) po _ UD ( q ) _ E ^ q) + ie 
x {2(P°-LU D (q)-ExM-P°j-^D(q)} , (48) 



where oooiq) = \/ q 2 + and E Sc (q) = y q 2 + m| c . The width of iV* g + — )■ J/^p is now 
given by 

r = Sm^ 1 '"*"* 1 (49) 
where tj/^ p ^ R means tj/^ p ^ R omitting the a ■ e operator. We take P° = M R = 4265 MeV 
and p = \ 1/2 (Ml,Mj^,M 2 )/2M R , while M p stands for the mass of the proton. By using 
the form factor of 48], we get 

r«_j/^ = 0.01 MeV , (50) 

with admitted uncertainties of the order of a factor two. Since T nN of the iV* g + (4265) was of 
the order of 3.8 MeV, now the cross section is about a factor 400 smaller than before. Yet, 
the fact that the background for J/ipp production is also smaller might compensate for it. 
But, from what we have said before, the cross section for r\ c p production is much bigger. 

On the other hand, for the resonances made out by VB, the J/ipp production cross 
sections are larger. One can repeat the calculations in this case. We sketch the derivation 
below. 
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J/lp 

p ,. M v p 




p , <r „ p 

FIG. 10: The pp — > J/ippp mechanism throughout the resonance N*^~ 

We shall make the estimate based upon the mechanism of the Feynman diagram of 
Fig. [TOj and we will consider the resonance N* 5 (AA18) coming from the interaction of vector 
mesons with baryons, one of which channels is J / ipp, which was considered in the Subsection 
III. C as seen in Table HXl By adding this new channel we found gxj/ipN — 0.85. Assuming 
the dominant decay channels of iV* as pA(For p°N, it should be added Cj = l/y/S) and 
dominance of the 7 term in the p°pp vertex, which goes then as (77 /a/2, and g = M v /2f, 
we obtain now 

rfgpp->;v e * 5 (44i8)p _ 9 2 M\ T pN Cj E(p')E(p) + p.p' + M 2 p' 

dcos9 4 s p° n (2M 2 - y/sE^p 1 ) + 2p.p' - M 2 ) 2 p 1 ' 

with p',p the p outgoing, incoming momenta in the center of mass frame, and p° n the 
p momentum in the iV*g(44i8) decay into pN. By means of Eq. (|5T|) and the width of 
iV* 5 (4418) — > J/ipp, we can calculate the cross section of the reaction pp — > J/ippp multi- 
plying the cross section of Eq. (!5T!) by the branching ratio of the resonance for the decay 
into J/ ipp. As one can see in Fig. [HI this cross section is of the order of 2 ~ 37 nb for a p 
beam of 15 GeV/c, depending on whether one includes or not the form factors. And for the 
dashed line, we also give the form factor for the NNp vertex and ^^(4418) as follows: 

A 2 - m 2 

p r p 

^ = 72 , , K 5 (53) 



^■N + (Pn* (4418) " "' 



2 

TV* (4418) / 



with A p = 1.3GeV and A^ = O.SGeV. 

This cross section is larger than the one we would obtain from the standard mechanism 
of Fig. [T21 which can be evaluated in analogy to the case of Fig. El Once again, using Eq. 
(EU) and r j/i^p of the resonance instead of Y^n we can obtain the differential cross section 
of the peak of the resonance: 6 ~ 50 nb/GeV 2 . 
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FIG. 11: The total cross section vs the p beam momentum for pp — > ppj/ip. The solid and 
dashed lines are calculated by p-meson exchange without and with form factors, respectively. The 
dot-dashed and dotted lines are calculated by Reggeon propagator with sq = 5 and 8 GeV~ 2 , 
respectively. 




FIG. 12: The standard pp — > J /ippp mechanism 



For the same reasons as for the iV* g (4265) production, we also consider the Reggeon 
exchange here. The Eq. ( 151ft becomes as follows: 

da pp ^ m8)p _ g*MlT pN Cl W)m + + j^ m M? t (54 ) 



where 



dcosO 4 s p° n v v ' y ' p 



a p (t)-l 



RJs,t) = at' (t)exp( — i— aJt)) t-t (55) 

py ' 2 » y> py 2 pW; cos( | ap(t) )r(^) + i) 1 ; 

a p (t) = 0.5 + 0.83t (56) 



a' p (t) = 0.83 (57) 
t = 2M 2 N - y/sE(p') + 2pp'cos9 (58) 
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When t — > mi, a a = 1 and R D ~ } , . For the estimation the value of Sn is taken to be from 

pi p p m A p —t u 

5 GeV -2 to 8 GeV~ 2 , such that \R p (s,t)\ 2 gives almost the same results as ^ m 2_ t ^2 when 
yfs < 3GeV. By using Reggeon propagator, the total cross section is about 0.008 ~ 0.06 nb 
for a p beam of 15 GeV/c. This is about two orders of magnitude smaller than the result 
by p-meson exchange. 

From the calculation above, we find that the cross section of this reaction is about two 
orders of magnitude smaller than that of the reaction pp — > pprj c , but it could be also 
appropriate to find iV*(44i8) because the J/ip has a large branching ratio to decay into 
lepton channels which are much easier to detect than hadron channels. 

Finally let us discuss the possibility of measurement of this reaction in the experiments. 
The PANDA (anti-Proton Annihilation at Darmstadt) Collaboration will study the pp reac- 
tion at FAIR, with the p beam energy in the range of 1.5 to 15 GeV/ c and luminosity of about 
10 31 cm _2 s _1 [52]. The range of the beam energy is very suitable to find the iV*(4265) and 
the iV*(4418), with cross sections estimated to be about 70nb and 2nb by the one-meson 
exchange propagators with off-shell form factors, which corresponds to an event produc- 
tion rate of 60000 and 1700 per day at PANDA/FAIR, or about lOnb and 0.02nb by the 
Reggeon propagators, which corresponds to an event production rate of 9000 and 20 per day 
at PANDA/FAIR. There is a An solid angle detector with good particle identification for 
charged particles and photons at PANDA/FAIR. For the pp — > pprj c reaction, if p and p are 
identified, then the rj c can be easily reconstructed from the missing mass spectrum against 
p and p. It is the same as the reaction pp — > pp.] ftp. So this reaction should be accessible 
at PANDA/FAIR. 

V. SUMMARY 

In summary, we find six states from PB and VB channels by using the local hidden 
gauge Lagrangian in combination with unitary techniques in coupled channels. All of these 
states have large cc components, so their masses are all larger than 4200MeV. The width 
of these states decaying to light meson and baryon channels without cc components are all 
very small. On the other hand, the cc meson - light baryon channels are also considered 
to contribute to the decay width to these states. Then rj c N and r] c A are added to the PB 
channels, while J / tpN and J / tpA are added in the VB channels. The decay widths to these 
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channels are not negligible, in spite of the small phase space for the decay, because the 
exchange D*(orD*) mesons were less off-shell than the corresponding one in the decay to 
light meson - light baryon channels. The total width of these states are still very small. We 
made some estimates of cross sections for production of these resonances at the upcoming 
FAIR facility. The cross section of the reaction pp — > ppq c and pp — > ppJ/ijj are about 
10 ~ 70nb and 0.02 ~ 2nb, in which the main contribution comes from the predicted 
A* g (4265) and A* g (4418) states, respectively. With this theoretical results, one can estimate 
about 9000 ~ 60000 and 20 ~ 1700 events per day at the PANDA/FAIR facility, respectively. 

The predicted N* 5 and A* E can be also looked for by many other processes, such as 
ep -> eN* £ at JLab's 12 GeV upgrade, Kp ->• A* E at JPARC, pp collisions, etc. 
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Appendix A: The C a \, coefficients 

In this Appendix we give the coefficients C a b in Eqs. f lTHj dU |2"5"|2"7|) for the several (/, S) 
sectors studied in this work. 

TABLE XI: Coefficients C ab in the Eq. (TJJ [55} for the PB system in the sector I = 3/2, 5 = 0. 









2 -1 1 
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TABLE XII: 
5 = 0. 



Coefficients C ab in the Eq. (JT31 E3 (27J) for the PB system in the sector I = 1/2, 





Z)S C DA+ 


rj c N 




TjN 


r/iV 






-1 




-1/2 


-1/V2 


1/2 


1 


DA+ 


1 


a/372 


-3/2 


1/V2 


-1/2 


1 



TABLE XIII: Coefficients C a6 in the Eq. (USES} for the PB system in the sector 1 = 1/2, S= -2. 









7rH 




77H r/'H 


KA 




1 


y/2 





V3/2 


1/V6 1/^3 


-V3/2 






1 





-3/2 


1/V2 1 


1/2 









^3/2 





-1/V3 l/x/6 






TABLE XIV: Coefficients C afe in the Eq. (HUES]) for the PB system in the sector / = 1, 5 = -1. 





A^c DE' C 




vrS 


7rA 




KiV K H 




v/2 











-1/V3 a/2/3 


-1 


DEI 


1 





1/V2 


-V3/2 


l/\/6 1/2a/3 


l/y/2 






1 


-a/372 


1/2 


-i/2 


a/3/2 
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TABLE XV: Coefficients C ab in the Eq. ([TBI 123] I2T |) for the PB system in the sector 1 = 0, 
S = -l. 





D S A+ 


DE C DE' C 


rj c A ttT, 


T]A 


j/A 




K H 


D S A+ 





-y/2 


1 


1/V3 


v/2/3 


-V3 





DE C 




-1 


1/V2 -3/2 


1/V6 


-l/2>/3 





a/3/2 


D~' c 




-1 


-y/372 ^3/2 


-1/V2 


1/2 





1/V2 
























TABLE XVI: Coefficients C ab in the Eq. CDS [25]) for the PB system in the sector I = 0, 5 = -3. 













2 



TABLE XVII: Coefficients C ab in the Eq. (ITU 125]) for the system in the sector I = 3/2, 5 = 0. 









L>*£ c pN K*T, 




2 -1 1 



TABLE XVIII: Coefficients C ab in the Eq. {HI [25]) for the VB system in the sector I = 1/2, 
5 = 0. 





D*Z C D*A+ 


pN 


ujN 


<f>N K*Z K*A 




-1 


-1/2 


V3/2 


10 


D*A+ 


1 


-3/2 


-V3/2 


1 
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TABLE XIX: Coefficients C ab in the Eq. CLUEIS for the VB system in the sector 1=1/2, S= -2. 





D*E' 

*^ s^c 




D*n c 


pE 


K*Y, 


loE 


0H 


K*A 


D*E' 

s c 


1 










V3/2 





-1/V2 


-V3/2 






1 








-3/2 







1/2 


D*Q C 









y/3/2 





y/S/2 









TABLE XX: Coefficients C ab in the Eq. (O [25]) for the VB system in the sector I = 1, S = -1. 









p£ 


pA 


wE 


iTW 


K*E 




\/2 














-1 -1 







1 





1/V2 


-V3/2 


-1/2 





1/V2 






1 


-a/372 


1/2 


\/3/2 





a/3/2 



TABLE XXI: Coefficients C ab in the Eq. flH [25]) for the VB system in the sector I = 0, S = — 1. 





L>*A+ 5*H C L>*H^ 


p£ ujA 


0A K*iV iT*~ 


5 S *A+ 


-a/2 





-1 -/3 




-1 


-3/2 -1/2 


a/3/2 




-1 


y/S/2 a/3/2 


l/>/2 



TABLE XXII: Coefficients C ab in the Eq. (HH [25]) for the VB system in the sector I = 0, 5 = -3. 













2 a/2 
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